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D ■ Abstract 



We present the correction to the decay rate into hadrons in the limit 
mtop » mz- 

> , 

00 1 Introduction 

O I Precision measurements of the decay rate into hadrons at LEP [|l| provide precise means 
^ to extract the QCD couphng constant from experiment. The approximation to the 

■ decay rate into hadrons is important for an accurate determination of a^. The hadronic 
^ ! decay rate is a sum of vector and axial-vector contributions of which the vector contribution 
Ph! is known to order from the calculation of ^^^(e+e^ j hadrons) 0]. The correctness 
of this calculation is strongly supported by 0] where the non-trivial connection between 
^ ■ the result and the approximation [H to deep inelastic sum rules was established. 

The axial-vector part of the hadronic Z^ decay rate was calculated to order in and 
confirmed in where renormalization group improvements were made. The Z^ decay 
^ I into 3 gluons in order has been calculated in 0. In this paper we present the total 
I correction to the Z^ decay rate into hadrons by calculating the axial-vector part in order 
al in the leading order of a large top mass expansion. 

2 Preliminaries 

For the Z^ decay rate into hadrons, the quantity to be determined is the squared matrix 
element summed over all final hadronic states. One can express this quantity as the 
imaginary part of a current correlator in the standard way 

J2<0\J^\h><h\r\0> = 2ImU^'\ (1) 
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W = ij d^ze'"-' < 0|T(J'^(z)J'^(0))|0 > = -g'^^'Iiiiq^) - q^^q'^U^iq^ 



(2) 



Here 



4cw ^ 



g=l 



(2/f) 



is the neutral weak quark current 



coupled to the bozon in the Lagrangian of the Standard Model. The hadronic decay 
width is expressed as 



^had — ^had + ^ 



had 



mz 



Imlli{m\ + 



(3) 



with the indicated decomposition into vector and axial- vector parts imposed by the struc- 
ture of the neutral current. 

Throughout this paper we use dimensional regularization P] and the standard modifi- 
cation of the minimal subtraction scheme the MS scheme [|lOl. For the treatment of 



the 75 matrix in dimensional regularization we use the technique described in |]TT| which 
is based on the original definition of 75 in We work in the approximation of 5 massless 
flavors and the top quark mass large compared to the Z^ mass. We should stress that the 



top quark does not decouple [|I2| from the axial-vector part due to diagrams of the axial 
anomaly type. 

The approximation for the vector part in effective QCD with 5 active massless 
quark flavors in the MS scheme reads 



q=l 



1 + — + 1.40923 — 

TT V TT 



a. 



QJ 

12.76706 ( — 

TT 



-.41318 ( — 

TT 



(4) 



with the Fermi constant Gp 



Here a. 



a 



(5), 



mz) is the coupling constant in 



effective QCD with 5 active flavors (the expression for af"* in NNL approximation can be 
found e.g. in 0). 

It is convenient to split the axial- vector contribution in a non-singlet and a singlet part 



had 



r 



A,NS 
had 



+ r 



A,S 
had ■ 



(5) 



The non-singlet part comes from Feynman diagrams where both axial vertices are located 
in one fermion loop. The non-singlet part can be reduced to the vector case by using the 
effective anticommutation property of the 75 matrix in the prescription that we use. The 
result in the effective theory with 5 active massless quark flavors is 



had ~ o /FT Z^V q / 

ov27r 



a. 



1 + — + 1.40923 { — { - 12.76706 (—^ ^ 

TT V TT / V TT 



(6) 
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3 The calculation 

The singlet contribution F^^^ comes from diagrams where each axial vertex is located in a 
separate fermion loop. The 3-loop and 4-loop singlet diagrams that we need to calculate 
are drawn in figure 1. 



Figure 1. The symbol is used to indicate an axial vertex. 



In the Standard Model quarks in a weak doublet couple with opposite sign to the Z° 
bozon in the axial part of the neutral current. That is why the contributions from light 
doublets add up to zero in the massless limit for singlet diagrams. The only non-zero 
contribution comes from the top-bottom doublet due to the large mass difference be- 
tween top and bottom quarks. The massless diagrams (i.e. without top quark loops) 
were calculated with techniques similar to the ones that were used for the calculation of 
o'tot{e~^e~ — 7 — hadrons) 0. The i?*-operation |jl3[ which subtracts both ultraviolet and 
infrared divergences was a necessary part of these techniques. For the calculation of mas- 
sive diagrams (i.e. with one or more top quark loops) we applied an asymptotic expansion 
in a large top mass m^. The theory of Euclidean asymptotic expansions was developed in 



14, 15 1 . For the large mass expansion of the diagrams we use the techniques developed in 



16| , p!7| . For the actual calculation we relied heavily on the symbolic manipulation program 



FORM [jl8[. Massless diagrams were calculated with the help of the package MINCER |19 |. 



Massive diagrams were calculated with a specially designed package for massive vacuum 
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3- loop diagrams which uses algorithms from ref. [Q. Both packages are based on the 
'integration by parts' method [^. Details of the calculation will be given in a longer 



paper where we will also present the higher order terms in the large mass expansion. The 
computations were done in an arbitrary covariant gauge. The cancellation of the gauge 
dependence in physical quantities gives a good check of the results. The result of our 
calculation in the leading approximation of the large rrit expansion for 6 flavors is 



.A,5 
had 



a 



(6)- 



a 



(6)^ 



TT 



71" 



(7) 



+CfTID (-1 + ICs - ll 



+ TlD 



108 ^ 18" ^ 12'''^\^) 6"'' 



Here rrit 



mt{fi) is the running top mass in the MS scheme. a^^\fi) = |^ is the coupling 



constant in QCD with 6 flavors, Cp = and Ca = 3 are the Casimir operators of the 
fundamental and adjoint representation of the color group SU{3), D = 8 is the dimension 
of the Lie algebra, Tp = ^ is the trace normalization of the fundamental representation, 
Nf = 6 is the number of quark flavors, ( is the Riemann zeta-function. Contributions with 



TC 



originate from terms containing In {—rriz — = {ln{'m 



ITT I 



when one takes the 

imaginary part of the correlator. Results for individual diagrams contain (2, C4 and Cs but 
these contributions add up to zero in the total result. Note that the order agrees with 
the known calculation [Q. 

The coefficients of the logarithms are in agreement with the required renormalization 
group invariance of the physical quantity 



d 
d^j? 



Of course the true physical quantity is T^ad- But from a theoretical point of view the 
singlet part Vf^f^ is renormalized independently of the non-singlet part and is therefore 
renormalization group invariant by itself. 

The approximation for the vector part and the axial non-singlet part are calculated 
in effective QCD with 5 active massless flavors where the top quark decouples. Although 
the top quark does not decouple from the singlet axial part we should convert it to an 
expression in the effective theory to give a consistent total Thad- The QCD decoupling 

~ Wi. The 



relations in two loop approximation in the MS scheme were calculated in 
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connection between the full coupling constant a^f^ and the coupling constant of effective 
QCD with 5 active flavors reads 



.(5) 



+ 



TT 



7r 



mt 



(9) 



Substitution of this expression in (|^ gives the expression for effective QCD 



+ 



TT 



5651 
" 216 



■^A,S 

had 



23 



8^/2' 



71 



31 



TT 



37 
12 

373 



+ ln{ 



rrir. 



+ 



mt 



+ Cs + + TZ^n{^) + — /n( 



36 



mt 



24 



12 



23, n,m 
-ln^( — 



mt 



12 



(10) 

where = is the coupling constant in effective QCD with 5 flavors. The term 

drastically diminishes the /U- dependence of eq.(|lO|) and makes it stable in a wide interval 
around /i k, mz- 

Summing equations (^), (|^) and (|10|) and putting the renormalization scale /u in the 
standard way equal to mz we get the hadronic decay width (in the approximation of 5 
massless quarks and a heavy top quark, in the leading order of the large mt expansion) 



had 



ov27r 



1 + — + 1.40923 

TT 



12.76706 



+ 



w) 



-.41318 



+ 



+ 



Gf^ y (2/(3))2 [l + ^ + 1.40923 (—)'- 12.76706 

V27T t^i [ V TT y 



a. 



TT 



+ 



37 
12 



+ 



31 
1^ 



,1Tir, 



-18.65440 + ^ln{'^) + ^In^ 



ml 



12 



m 



(11) 

Here and below mt = mt{mz) is the MS top mass at the scale mz- One may relate it to the 



pole mass through the expression mt{mz) 



m. 



pole 



ln(- 

V 1 



pole 



+ 



0(al 



which is known in the NNL approximation or relate it to mt{mt) through the expression 



mtlmz 



mtlmt 



as{mz) 



rnf{mt) ■ 



+ 0{al 



. This would correspondingly modify the 



coefficients of the a? term. 
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Substitution of the values of the physical parameters into eq.([Tl|) from gives (in 
the MS scheme) 



THad{GeV) = 1.671 



1 + 



+ 



.9502 + .1489 /n(^)) + 
mi I 



as 

+ { — 

n 



rrir 



-15.650 + .2564 ln{-^) + .2853 ln\-^ 



mi 



mi 



'12) 



Note that the logarithms in the term tend to cancel each other for mt{mz) ~ 140Gey. 
For this value of mt the calculated order a^ singlet contribution to Thad adds about 20% 
to the previously known sum of vector and axial non-singlet contributions of order a^. 

From eq.(0) one can also obtain the approximation to o"to((e+e^ 1-,Z^ —>■ 
hadrons) in the energy range below the top quark threshold. 
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